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Abstract 

Given two linearly independent matrices in so(3), Z\ and Z2, every rotation matrix 
Xj G S'0(3) can be written as the product of alternate elements from the one dimen- 
sional subgroups corresponding to Z\ and Z2, namely = e^^^^e^^'^e^^^^ • • • e^^*". 
The parameters ti, i = 1, s are called generalized Euler angles. 

In this paper, we evaluate the minimum number of factors required for the factor- 
ization of Xf G S0{3), as a function of Xj, and provide an algorithm to determine 
the generalized Euler angles explicitly. The results can be applied to the bang bang 
control with minimum number of switches of some classical control systems and of two 
level quantum systems. 

Keywords: Decompositions of Lie groups, Rigid Body Dynamics, Geometric Control, Two 
Level Quantum Systems. 

1 Introduction 

In this paper, we deal with the problem of steering control for bilinear systems of the form 

X = Ax + Bxu, (1) 

where x G R m is a control function, and A and B are skew-symmetric 3x3 matrices, 
namely matrices in so(3). Several systems in applications have the structure (|I]). In par- 
ticular the most common example is given by the dynamics of the rigid body where one 
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component of the angular velocity is seen as the control u and the others are held con- 
stant. The fundamental matrix of equation (|T|) represents the orientation of the rigid body. 
Another example is the lossless electrical network dealt with in 0. A two level quantum 
system driven by a single time varying component of an electro-magnetic field also has the 
structure (Q) where x represents the state G (E ^ and the matrices A and B are in the 
Lie algebra su{2). Because of the connection between the Lie groups SO{Z) and SU{2) the 
results presented here can be applied to the latter system as well. 
The fundamental matrix of the system (P, X, satisfies 

X = AX + BXu, (2) 

with initial condition equal to the 3x3 identity matrix. It follows from the results of § 
that if A and B are linearly independent (and therefore generate so(3) which has dimension 
3 and no two-dimensional subalgebras) a piecewise constant control is sufficient to steer the 
state of @j from the identity to every matrix Xf in 5*0(3) and, as a consequence, the state x 
of (|l]) between two states with equal length. Let us assume now that the control u is allowed 
to attain only two values, M and N . Define Zi := A + BM and Z2'-= A + BN, and assume 
a factorization of the desired target state Xf of the type 

is known with ti,t2,...,ts > 0[|. Then a piecewise constant control equal to M for time 
tg, N for time t^-i, M for time t<j_2 and so on, drives the state of (0) from the identity 
to Xf in (^). This idea, involving Lie group decompositions, has recently been used to 
prescribe controls for quantum mechanical systems where the underlying Lie group is the 
group of special unitary matrices of dimension n, SU{n) (see e.g. 0, ||12|, [l^, and 



references therein). If the control is bounded in magnitude, namely \u\ < M we can choose 
N := —M and M as the two values for the control. From a practical point of view one would 
like to have a factorization of X/ in terms of the matrices Zi and Z2 that involves the least 
number of factors, so that the control law has the minimum number of switches. Moreover, 
an algorithm is needed to evaluate the generalized Euler angles tj, i = 1, s. This paper 
is devoted to the solution of these two problems. Constructive factorizations of SU (2) and 
5*0(3) can be found in the papers [|T3] that, however, do not consider the problem of 
minimizing the number of factors. 

The paper is organized as follows. In the next section we give some preliminary defi- 
nitions that will be used in the following and recall some results proved in |10| concerning 
factorizations of elements of the Lie group 5*0(3) of the type (|]). We also transform, using 
a change of coordinates, every pair of linearly independent matrices G so{3) into a canonical 
form, that will be used in the following sections, without loss of generality. In Section 3 we 
evaluate the minimum number of factors needed in a factorization of matrix Xf G 5*0(3) of 
the type (^ given Zi and Z2. In Section 4 we give an algorithm for the determination of 

^This is done without loss of generality since the one parameter subgroups corresponding to Zi and Z2 
are closed, namely the functions e^^'^* are periodic. 
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the generalized Euler angles ti, ...,ts- We discuss applications to the control of classical and 
quantum systems in Section 5. 



2 Preliminaries 

The inner product < •, • > between two elements of so(3), Zi and Z2 is defined as 

< Zi, Z2 >= Trace{ZiZ^). (4) 

If < Zi, Z2 >= 0, the maximum number of factors s needed to express a matrix Xf in 5*0(3) 
as in (0) (maximum over 5*0(3)) is three, and the factorization in (^ is the classical Euler 
resolution of a rotation (see e.g. [|l4l) (modulo a change of coordinates and a re-scaling of 
the variables t). The parameters ti are called Euler angles and their calculation is standard 
matter (see e.g. 11^, pg. 297). 



In |T^, it was shown that, for every pair of matrices Zi, Z2, the number of factors needed 
to express an element Xf G 5*0(3) is uniformly bounded, over 50(3) (see also [@] and ||15| 



for generalizations to every compact Lie group). The maximum value for s (maximum over 
5*0(3)) is called the order of generation of 50(3) with respect to Zi and Z2. It has been 



calculated in [rU| and it only depends on the value of the cosine of the angle between Zi and 
Z2, namely 

Sll^ ,5) 

< Zl, Zl > 2< Z2, Z2 >2 

If ip = 0, the order of generation is equal to 3 and we obtain the classical Euler resolution 



of a rotation. Our treatment in the following was inspired by the proof in However 



most of the treatment in |10| is carried out using stereographic projections and translating 
the problem to the induced subgroup of the Moebius group. We shall treat the factorization 
of every element in 5*0(3) by working on the unit sphere in R ^ and looking at 5*0(3) as a 



transformation group on the sphere [11 



We now show that there is no loss of generality in assuming that Zi and Z2 in (^ have 
a special form which we shall describe. We shall call Shk, h < k, the matrix in so(3) which 
has zeros everywhere except in the h, k-th [k, h-th) entry which is equal to 1 (—1). Given a 
matrix Zi, there exists a matrix Ti G 5*0(3) such that 

TiZiT[ = Ai5i2, (6) 

Ai 7^ 0. This can be easily seen by choosing Ti := [t>i, t>2, fs]"^, with ^3 such that vJZi = 
and with norm equal to one and Vi and V2 such that {vi,V2,Vs} form an orthonormal basis 
in R ^. We also set 

TiZ2Tf := aSi2 + 65i3 + c523. (7) 

Choose now T2 := c^^^^ with 6 such that hcos{6) + csiniO) = 0, with b and c given in (0). 
Then we have 

T2TiZiT^T2 = Ai5*i2, (8) 
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T^TrZ^T^T^ = aSu + ^^23, (9) 

for some parameter d ^ 0. Therefore, we can always assume that, in appropriate coordinates, 
the matrices Zi and Z2 have the form Zi := XiSu and Z2 := aS'12 + dS23, respectively. 
Moreover we can divide Zi by Ai 7^ (this has the only effect that, in the matrices of the 
form e^^*, t has to be scaled by a factor Ai) and analogously we can divide Z2 (in the new 
coordinates in (^) by d 7^ and therefore the parameter t in the subgroup e^^* has to be 
scaled by a factor d. Define p := ^. We can assume, without loss of generality, that the 
matrices Zi and Z2 are given by 

Zi := Su, (10) 

and 

Z2 := pSi2 + S23, (11) 

and we shall do so in the following. Notice that the above manipulations do not modify the 
value of the parameter ip in (|^) which is given, in terms of p, by 

^ = ^=^- (12) 

3 Decompositions with minimum number of factors 

Assume now that an element Xf G 5*0(3) is given, to be expressed as in (|]), with Zi and Z2 
given in (|l3), (|ll|)- We give in this section a procedure to determine the minimum number 
of factors necessary as a function of Xj. 

We assume p in (|llD different from zero (the case p = corresponds, from ([T^), to 



Zi and Z2 orthogonal to each-other and therefore the decomposition is the standard Euler 
decomposition). Define two sequences {zk} and {fk}hjzo = fo = —l 

fk:=^^^l-zl + Zk (13) 

Zk+i ■= — — j/fc - Zk. (14) 
1 + p'^ 

We have the following Lemma. 

Lemma 3.1 There exists an index k >1 such that fk is defined (\zk\ ^ I) for every k < k, 
fk < 1, for every k <k and > 1. 

Proof. First notice that if |p| < 1, the Lemma is true with k = 1 since < 21 < 1 and 
/i > 1. Let us assume |p| > 1. We first show that fk well defined and /a; < 1 implies that 
fk+i is well defined, namely that |2;a;+i| < 1- Then we show that there exists the first value 
of k, k, such that > 1. 

Assume fk < 1. From (|l^), we obtain 

^/i^k<m-zk), (15) 
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which gives, taking into account \zk\ < 1, 



1 + P 

Consider the expression of Zk+i obtained combining (p!3|) and (|14| 

1 + 



l<Zk< ^. (16) 



Zk+i = — — -{ — ^l-zl + Zk)-Zk. (17) 



Consider function of in the interval defined in ([iq) . This function is always 

increasing from the value z^+i = at 2;/,. = —1 to the value Zk+i = 1 at z^ = ^p^- In 
particular we always have < 1 which implies that fk+i is well defined. To show the 

existence of ak such that > 1, we evaluate Zk+i — Zk using ([13|) and (0). We obtain 



Zfc+i -Zfc = Y:^i\p\\/'^- 4- Zk)- 



Using the second inequality in (p!6|), we obtain 



2\p 

and plugging this into (|TB]), we obtain 



1 - > (19) 



Zk+1 - Zk> ^ , , - Zfe) > :;— — ^, (20) 



where, in the last inequality, we used inequality (|T^) again. Therefore the sequence {z^} is 
increasing by at least at each step and since fk > Zk for every A;, from (p!3|), we must 
have a value of the index k such that fk> 1. This concludes the proof of the Lemma. □ 



We now relate the finite sequences {zk} and {fk}, k = 0,1, ...,k defined in (0) and (|14D 
to the minimum number of factors needed in the factorization (^). Consider a given target 
matrix Xf := {xij}, i,j = 1,2,3 to be factorized. We define a function 0{Xf) which is 
equal to 1 if ^3^3 = —zq = 1, it is equal to 2 if < "2:3^3 < Zi and Xi^3 = p(— a;3^3 + 1) and 
equal to 3 if < —^3,3 < and Xi^s 7^ p(— 3^3,3 + 1)- In cases not considered above, let k 
be the highest value of the index k such that 

Z-k < -3^3,3 (21) 



(recall from (20) that Zk is increasing at each step by at least a given amount). Then we 
have 

0{Xj) = 2^ + 2 if sign{p)x,;, > -\p\{xs,3 + ft) (22) 

OiXf) = 2k + 3ii stgn{p)xi,3 < -|p|(x3,3 + fk). (23) 
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The following Lemma gives the minimum number of factors in the factorization @ assuming 
that the first factor on the right is of the form e^^*'' with tg > 0. The proof of the Lemma 



reveals the geometric meaning of the finite sequences {zk} and {fk} defined in (|T3D, (|Tj). 
We denote the minimum number of factors needed to express a general matrix as in 
by MIN{Xf). 

Lemma 3.2 Assume that Xf is such that the factorization with minimum number of factors 
in (j^ starts with a nontrivial factor of the type e^^* on the right. Then 

MINiXf) = OiXj). (24) 



Before giving the proof of the Lemma, we describe the geometry of the above construction. 
Considered as a transformation on the sphere of radius 1 centered at the origin, Xf transforms 
the South Pole Pg := [0,0,-1]-^ into a point Pj := [— ^1^3, — 0:2,3, —Xs^s] (which is just the 
negative of the third column of Xf). Conversely, any matrix Xf such that Pf = XfPg is 
equal to Xf up to a factor that leaves Ps unchanged. Such factor will in general have the 
form e^^* (recall (|l3)) and therefore we have Xf = XfC^^^. We would like to find any product 
with minimum number of factors 

Xf ■= e^i*ie^2t2g^ii3 . . . gZ2<.-ig^it.^ (^25) 

(with ts possibly equal to zero) such that Pf = XfPg and then to obtain Xf as Xf = 
XfC^^^" . From the assumption that the minimum number of factors for Xf is obtained with 
a nontrivial factor e'^^*" on the right, the minimum number of factors will be given by s. 
This observation can be interpreted in the language of coset spaces and homogeneous spaces 
(see e.g. |jll[). The subgroup H := {X G SO{?>)\X = e^^*,t G R } is the isotropy group of 



the South Pole Pg, namely the set of all the elements of 5*0(3) that leave Ps fixed. There 
exists an isomorphism between elements of the coset space S0{3)/H and elements of the 
sphere S^. In the expression (|^) we use the last term e^^*'' to move inside a coset while the 
remaining factors are used to go from one coset to the other, namely from one point on the 
sphere to the other. We now look for a transformation Xf in (pSf) transforming Pg to Pf 
with minimum number of factors. 

On the sphere S^, every element of the form e^^* corresponds to a rotation about the z 
axis. Each point on the sphere 5*^ follows a trajectory on a circle which is the intersection 
of the sphere with a horizontal plane. The value of the z coordinate of the point is not 
changed by this rotation. Every matrix of the form e^^* corresponds to a rotation about 
the axis defined by the vector Up := [1,0, p]^ (points on the line through the origin parallel 
to this vector are left invariant by the rotation). Under the action of this rotation, every 
point on the sphere follows a trajectory on a circle which is the intersection of a plane 
perpendicular to ftp and the sphere 5''^. Each such plane forms an angle 9 := tan~^^ with 
the X — y plane. If we consider a trajectory e^'^^Ps := [x{t).,y{t)., z{t)Y ■, the maximum value 
for the coordinate z{t) will be obtained at zi defined in ([I^ ) (when t = it). Let us call this 
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point Pi. Following a horizontal trajectory e^^*Pi, for t = tt, we obtain a point which is 
opposite to Pi- Let us denote this point by Qi. Following from Qi a trajectory e^^^Qi again 
up to e^^'^Ql, we obtain a point with z coordinate given by Z2 in ([T^). The value /i is the 



2;— coordinate of the intersection of the plane perpendicular to ftp, containing the point Qi, 
and the 2;-axis. Notice that Z2 is the maximum value that can be obtained for z starting 
from Pi and with just one switch from one type of trajectory to the other. Continuing this 
way one obtains the elements of the sequences {z^}, {fk}- It follows from this geometric 
description that {zk} is an increasing sequence and it was proven in Lemma 3.1 that it is a 
finite sequence (See also the Remark following the proof of Lemma 3.2). Figure 1 describes 
(in a two-dimensional plot) the trajectories on the sphere. In this Figure, k, defined in 
Lemma 3.1, is equal to 5. We have denoted by P^, k = 0,1, ...,5 the points on the sphere 
with z— coordinate equal to z^, k = 0,1, ...,5. denotes the point whose 2;— coordinate if 
fk, k = 1, 5. We now use this Figure to complete the proof of the Lemma. 

Proof of Lemma 3.2. We shall refer to Figure 1 and the above discussion. Let Pj := 

[—2^1,3? —^2,3, —3^3,3]"^- If 2^3,3 = 1 = —Zq, then Xf is of the type e^^^ and clearly MIN{Xf) = 
1. If zq < —Xs^s < Zi and Xi^s = p(— ^3^3 + 1) then Pf belongs to the intersection of the 
plane x + p{z + 1) = with the sphere S^. The point Pf can be reached by (possibly) 
following a trajectory of the type e^^^ (which leaves P^ := [0,0,-1] unchanged) followed 
by a trajectory of the type e^^*. In this case, since we have assumed that the last factor 
on the right in (|^) is a nontrivial e^^* factor, we have MIN{Xf ) = 2. Analogously, it is 
easily seen that MIN{Xf) = 3 if < —0:3^3 < zi and Xi^3 7^ p(l — ^3^3). Now notice 
that to reach a point with z coordinate z, with Zk < z < Zk+i we need to cross the circle 
Ck := {{x,y, z)\z = Zk,x'^ + y'^ + z^ = 1, {x, y, z) 7^ (x^, yk, Zk)}. In order to cross any point of 
Cki the minimum number of pieces of trajectory (including possibly the first one of the type 
e^^*, if assumed nontrivial, and the last one to leave Ck)-, is 2k + 2. This is clear when k = 1 
and follow by induction for the other values of k, noticing that we must cross Ck-i before 
crossing Ck- To reach a point Pf := [— 2:1,3, — 0:2,3, —3:3,3]"^ such that zj, < — ^3,3, with k <k, 
we need to cross Cj^ and the minimum number of factors to do that is 2k + 2. No other factor 
is needed if Pf = [—3:1,3, —3:2,3, —3:3,3]-^ is below the plane with equation x + p{z — Z^) = 
while another factor is needed if Pf is above this plane. This accounts for the inequalities 

□ 



Remark: It is possible to show that the sequence {zk} in (|I^ ([I^ can be obtained by 
Zk = —cos{k[3), for some angle /3 obtained as /3 = cos~^zi. From a geometric point of 
view, j3 is the angle in the y — z plane between the segments ODk and ODk+i, where O 
denotes the origin and Dk : e^^^Pfc, k = 0, 1, 2, 5. The points Dk are the midpoints of the 
lines representing a circle Ck in Figure 1. This angle is the same for every k. This gives a 
geometric interpretation and an alternative proof of Lemma 3.1. □ 

The above Lemma solves the problem of finding the minimum number of factors to 
express Xf in the form (0) if we assume that the last term on the right is of the type e^^*. 
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Figure 1: Geometric Construction for Lemmas 3.1 and 3.2 
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This assumption can be relaxed by considering a change of coordinates T, 

zR^ 1 



T : = 





1 



1 







p 



We have 



TZiT^ 



1 + 
1 

■■Zo. 



(26) 

(27) 
(28) 

(29) 

with a term of the type e^^^ first on the right, is the optimal factorization. Then a factoriza- 
tion with a term e^^* first on the right and s factors is the optimal factorization for TXfT^ 
and viceversa. Therefore according to Lemma 3.2 we have MIN{Xf) = 0(TXfT'^). We 
conclude with the following Theorem. 



Assume that 



Xf 



Theorem 3.3 



MIN{Xf) = min{C(X/), 0{fXff^)}, 



(30) 



with T given in (\2E). 



4 Evaluation of the Generalized Euler Angles 

The geometric analysis of the previous section gives a method to determine the generalized 
Euler parameters corresponding to the optimal factorization. Let us assume, without loss of 
generality, that MIN{Xf) = 0{Xf) namely, the optimal factorization has a nontrivial term 
of the type e^^* last on the right. Referring to Figure 2, we have labeled each region with 
a number denoting the minimum number of factors needed to drive Pg to Pf in that region 
(Including the last nontrivial factor on the right of the type e^^*). If Pf := — [xi^s, X2,3, Xs^s]"^ 
is in an odd region, such as Po in Figure 2, (namely it is strictly above a plane dividing 
a region between two planes z = constant, except for the Region 3, which includes points 
strictly below the plane x + p{z + 1) = as well) then an optimal factorization for Xf is 

Xf = e^i*ie^2*2e^^^e^27r . . . ^Zm ^Z27r ^z^ts _ ^^^^ 

We first determine so that, defined L := e^^^^e^^^'e^'"'^ ■■■e^^'^e^^'' := {kj}, = Xg^g. Then 
we determine ti so that e^^^^LPg = Pf, where Pg denotes the South Pole Pg = [0,0,-1]"^ 
and then tg such that e^^*^Le^^*° = Xf. Notice that each step involves the evaluation of just 
one parameter. Notice also that the optimal factorization is not unique and, in the above 
factorization, we could have, for example, replaced the term e^2*2g^i7r -with a term e'^2*2g^iti 
for appropriate values ii and t2 (see the alternative path with bold face lines in Figure 2). 
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Figure 2: Optimal paths on the sphere 

If MIN{Xf) is even (Py = Pe hi Figure 2) then we have that the optimal factorization is 
given by 

Let, in the sequence (0), QT^), Zk be the largest value of Zk such that Zk < —3:3,3. Then, 
we consider a point P := [x,y,z]'^ intersection of the planes (x + xi,3) + p{z + 0:3,3) = 0, 
z = Zk and the sphere + + = 1. Then we determine t2 so that, defined L : = 
gZitagZaTTgZiTT. ..gZiTTgZaTrp^^ j^^ve LP^ = P. Then we determine ti so that e^^tip = and 

finally we determine tg so that e^^^^Le"^^** = Xj. In this case too, the optimal factorization 
is not unique. 

5 Applications 

The results of this paper can be used to prescribe bang bang type of controls for bilinear 
systems whose state varies on the Lie group S0{3), with minimum number of switches. 
This technique of control can be applied to the dynamics of a rigid body where the angular 
velocity is seen as control. The same technique can also be employed for the control of 
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switched electrical networks 0, with minimum number of switches. 

In recent years there has been a large amount of interest in the control of systems of the 
form (|l|) with A and B in the Lie algebra su{2). This Lie algebra is isomorphic to so{3). 
These systems model the dynamics of two level quantum systems with just one control 0. 
Constructive factorizations of elements of the Lie group SU{2) have been given in 0, 



131 , and used for control. In particular the factorization of p] gives a worst case number of 
factors which is greater than the minimum by at most one. The algorithm presented in this 
paper can be used to determine the optimal factorization for elements Xj in SU{2) as well 
and therefore to prescribe a control algorithm for two level quantum systems with minimum 
number of switches. 

Let (f) denote the isomorphism between su{2) and so{3) which maps the Pauli matrices 

to 25*1^3, 252,3, —25*1,2 respectively. Let Zi and Z2 be two linearly independent matrices in 
su{2). We look for the factorization of Xf of the type 

Xf = e^'^'e^'^^ ■ ■ - 6^'^% (34) 

with minimum number of factors. The isomorphism between su{2) and so(3) induces a 
homomorphism between elements of the corresponding Lie groups, : SU{2) — > 50(3), 
which is given, if 5 = G 5f/(2), by 0(5) := e''^*^^\ This homomorphism is two to one in 
that to ±5 in SU{2) corresponds the same element in 50(3) (for a more detailed treatment 
of the relation between the Lie groups SU{2) and 50(3) see e.g. IjT^. See also for 
applications to control). Let Xf be the element in 50(3) corresponding to Xj under this 
homomorphism and Zi and Z2 the elements of so(3) corresponding to Zi and Z2. If 

X/ = e^i*ie^2*2 ■ ■ -6^^*% (35) 



is the optimal factorization for Xf then s is the optimal number of factors for Xf in ( p4D 
as well. The generalized Euler parameters can also be easily determined. If we use the 
same values for ti,...,ts in (|35|) and (|3^) we obtain a matrix which is ± the desired Xf. 



This affects the quantum mechanical state x in (|T]) by an overall phase factor which has no 
physical meaning. In any case, the minus sign can be easily eliminated by changing the value 
of just one of the parameters so as to change one factor 5 into —5. This is always possible 
since each one dimensional subgroup in SU{2) that contains 5 also contains —5. Therefore, 
we can find an optimal factorization for any element in SU{2) as well. This can be easily 
extended to any Lie algebra isomorphic to su{2) and the corresponding Lie group, which is 
known to be isomorphic to either 50(3) or SU{2) 0. 
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